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We present searches for the leptonic decays B^ t^v and the Icpton flavor violating decays 
i?" l^T^ , where I = e,n, with data collected by the BABAR experiment at SLAC. This search 
demonstrates a novel technique in which we fully reconstruct the accompanying B in T{AS) — > BB 
events, and look for a monoenergetic lepton from the signal B decay. The signal yield is extracted 
from a fit to the signal lepton candidate momentum distribution in the signal B rest frame. Using 
a data sample of approximately 378 million BB pairs (342 fb~^), we find no evidence of signal 
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in any of the decay modes. Branching fraction upper limits of B{B^ e^v) < 5.2 x 10 ^, 
B{B+ ^+v) < 5.6 X 10-^ e+r") < 2.8 x lO"''^ and ^ m'*'^"") < 2.2 x 10"^ are 

obtained at 90% confidence level. 

PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 



In this paper, we present searches for the decays 
t^v and the lepton flavor violating decays — > ^^t^, 
where £ = e,^ [H, using a technique in which the ac- 
companying B in the event is exclusively reconstructed. 
This method has not previously been used for searches 
for these modes and, although statistically limited with 
the present BaBAR data sample, shows promise for future 
studies at, for example, a high luminosity Super B fac- 
tory 0. While the former decay modes are allowed in 
the Standard Model (SM) and the latter are not, both 
are potentially sensitive to New Physics (NP) effects, 
such as contributions by neutral and charged non-SM 
Higgs 

Searches for rare B decays with neutrinos in the final 
state are challenging due to the limited availability of 
kinematic constraints. However, purely leptonic B de- 
cays involving an electron or a muon have a clear ex- 
perimental signature in the form of a high momentum 
lepton. Combined with clean theoretical predictions due 
to the lack of QCD contributions in the final state, such 
leptonic B decays present an ideal place to test the SM 
against NP models. 

In the SM, B^ l^v decays proceed via an annihila- 
tion of b and u quarks into a virtual boson. In the 
SM the branching fraction for this type of decay is given 
by i: 



B{B+ tvt) 



"^1 1 f 



2 IT/ 

Vub\ TB 



(1) 

where Gp is the Fermi coupling constant, m; is the 
lepton mass and m^, tb and /s are the mass, life- 
time and decay constant for the B meson. \Vub\ is the 
Cabibbo-Kobayashi-Maskawa matrix element which de- 
scribes the transition from 6 to u quarks Q . Within the 
SM, a determination of any one of the leptonic branch- 
ing fractions represents a determination of the product 
\Vub\'fB, which can be directly compared with determina- 
tions from lattice calculations f6| , _B-mixing and semilep- 
tonic decay measurements [7, 8]. As seen in Eq.(IT]), the 
decay rates are proportional to m/^, resulting in SM pre- 
dictions for the /i and e modes which are suppressed by 
factors on the order of 250 and 10^, respectively, com- 
pared with the r mode. Taking the branching fraction 
B{B+ T+Vr) = (1.31 ±0.48) X 10"'' from the combina- 
tion of recent BaBAR and BELLE results 0, [l3| implies 
Bsm{B+ ^ M+i/^) - 5.2 X 10-7 and Bsm{B+ e+i^e) ~ 
1.2 X IQ-^^. New Physics contributions to these pro- 
cesses can enhance or suppress the decay rates compared 
to the SM, and may either preserve or violate the rela- 



tive rates of the three leptonic modes depending on the 
particular NP model [1, 11|. Thus, the e and /j, modes be- 
come particularily interesting in light of recent evidence 
for the -B^ — > r+ Vr decay mode. Currently, the most 
stringent published limits on B^ i^v are from the 
BELLE collaboration with e+v) < 9.8 x 10"'^ 

and B{B+ ^i+v) < 1.7 x 10"*^ 12]. Earlier stud- 
ies by CLEO and BaBAR collaborations are also avail- 
able In, El. 

The lepton- flavor-violating (LEV) leptonic B decays, 
such as B° £^ , are forbidden in the SM in the ab- 
sence of non-zero neutrino masses, but can occur via one- 
loop diagrams if neutrino oscillations are included. The 
rates of such processes, however, would be substantially 
below current or anticipated future experimental sensi- 
tivities. On the other hand, many models of physics be- 
yond the SM, in particular supersymmetric seesaw mod- 
els d, predict dramatically higher rates for these de- 
cays. In the case of Higgs-mediated LEV processes, cou- 
plings to heavier leptons are favored, making B'^ £^ 
particularily interesting. In the general flavor-universal 
MSSM, the branching fractions allowed for 5° £^t^ 
are ~ 2 x IQ-^" Such decays could be within the 
reach of a Super B factory with a data sample of 10 to 
50 ab"^. The current best experimental limits on the 
branching fractions for these two decays are B{B'^ 
e+r") < 1.1 X 10-4 and B{B" -> ^i+t-) < 3.8 x 10-^ 
set by the CLEO collaboration with 10 fb-^ of data [l^. 

The searches described in this work are based on a 
data sample of approximately 378 million BB pairs, cor- 
responding to an integrated luminosity of 342 fb-^ col- 
lected at the T(45) resonance by the BaBAR detector 
at the PEP-II asymmetric e^e~ storage ring. Recon- 
structing the accompanying B meson in specific hadronic 
modes prior to the signal selection allows the missing mo- 
mentum vector of the neutrino(s) to be fully determined. 
The resulting increase in the energy resolution and the 
ability to infer the signal B meson rest frame provide the 
extra kinematic handles that permit signal events to be 
cleanly distinguished from the background. Previous B 
factory searches for 5+ £'^v and B^ have 
used an inclusive method in which the accompanying B 
is not explicitly reconstructed. This results in a signifi- 
cantly higher efficiency, but also a substantially increased 
background compared with the exclusive reconstruction 
method presented here. With the current level of lu- 
minosity, the inclusive method provides more stringent 
limits. However, due to the very low background achiev- 
able with the exclusive method, the two methods have 
about the same sensitivity for a statistically significant 
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observation. The method described in this paper will be 
the preferred approach for the high-precision studies of 
leptonic B decays. 

Charged-particle tracking and dE/dx measurements 
for particle identification are provided by a five-layer 
double-sided silicon vertex tracker and a 40-layer drift 
chamber contained within the magnetic field of a 1.5 T 
superconducting solenoid. A ring-imaging Cherenkov de- 
tector provides efficient particle identification. The en- 
ergies of neutral particles are measured with an electro- 
magnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals arrayed in a cylindrical barrel and in a forward 
endcap. Muon identification is provided by resistive plate 
chambers (partially replaced by limited streamer tubes 
for a subset of the data that is used in this analysis) inter- 
leaved with the passive material comprising the solenoid 
magnetic flux return. Signal efficiencies and background 
rates are estimated using a Monte Carlo (MC ) si mula- 
tion of the BaBAR detector based on GEANT4ll3|- The 
BaRAR detector is described in detail in Ref. [l7| . 

Reconstructed charged tracks are assigned a particle 
hypothesis based on information from detector subsys- 
tems, ifg candidates are selected by combining op- 
positely charged tt candidates and requiring that the 
TT+TT^ invariant mass satisfies 0.47GeV/c^ < m^r+Tr- < 
0.52GeV/c^. 7r° candidates are obtained from the com- 
bination of EMC clusters with no associated tracks, each 
with a T{AS) center-of-mass (CM) rest frame energy 
greater than 20 MeV, for which the 77 invariant mass 
satisfies 115 MeV/ < m^^ < 150 MeV/ c^. 

Over 96% of the tirne, the T{AS) resonance decays into 
a pair of B mesons [18]. Since the CM energy is pre- 
cisely known at PEP-II, exclusive reconstruction of one 
of the two B mesons, which we denote -Btag, fully de- 
termines the momentum four- vector of the other B me- 
son in the event. Charged and neutral B meson can- 
didates are reconstructed in hadronic final states of the 
form B D^*''Xhad- The reconstruction procedure be- 
gins with a Z3(*)o or seed, to which charged and 
neutral pions and kaons (which form the Xhad system) 
are then added. The combination of the D^*-' and Xhad 
with the lowest value of AE = \Eb — -Bbeaml that sat- 
isfies the condition AE < 0.2 GeV is chosen as the -Btag 
candidate, where Eb is the energy of the reconstructed 
B meson and ii'beam is the beam energy, both evaluated 
in the CM frame. We reconstruct D*^ in the £)+ tt^ and 
D° 7r+ channels, and D*° in the D° tt" and D° 7 chan- 
nels. The is reconstructed in the modes K~ 7r+ 7r+, 
K'^n+, K^TT+ir", R- 7r+ tt+ tt" and K'^n+n+ir- . For D" 
we consider the modes 7r+, 7r+ tt*^, 7r+ 7r+ 
w~ and if^TT+TT". 

Although multiple Z)'-*-'A"had combinations may be 
present in a single event, this procedure permits, at most, 
a single i?tag candidate to be retained in any given event. 
For the Stag candidate, we define the energy substituted 
mass. Toes = V -^bcam ~ Pb > where pB is the momentum 



of the i?tag candidate in the CM frame. i?tag candidates 
that are correctly reconstructed peak in toes near the 
nominal B meson mass, while incorrectly reconstructed 
-Btag candidates produce a combinatorial distribution. 
The signal events are required to lie within the range 
5.270 GeV/c^ < toes < 5.288 GeV/c^. This reconstruc- 
tion procedure results in a yield of approximately 2500 
(2000) correctly reconstructed B^ (-B°) candidates per 
fb~^ of data. 

Because the two B mesons are produced with very 
little momentum in the CM frame, BB events typi- 
cally produce a more isotropic distribution of particles 
in the detector than non-resonant ("continuum") back- 
grounds. Such backgrounds {e~^e~ ff, where / repre- 
sents u,d,s,c or any charged lepton) are suppressed by 
requiring R2 < 0.5, where R2 is the ratio of the second to 
the zeroth Fox- Wolfram moment [l^l computed using all 
charged and neutral particles in the event. Further sup- 
pression is achieved by requiring |cos0t| < 0.90, where 
6t is the angle between two thrust axes in the CM frame, 
the first computed using the particles from the -Btag , and 
the second using all other particles in the event. 

All particles that are not used in the -Btag reconstruc- 
tion are considered candidates to be included in the re- 
construction of the signal B meson. Since the CM en- 
ergy is precisely known, reconstruction of the -Btag fully 
determines the -Bgignai 4-vector. This permits the 2-body 
kinematics of the signal decays to be exploited. In partic- 
ular, these decays are expected to contain an electron or 
a muon with a momentum p*, in the -Bgignai rest frame, 
of about 2.64GeV/c (2.34GeV/c) for the B+ £+v 
(-B" — > i^T") channels, very close to the kinematic end- 
point for B decays. 

Signal candidate events are initially selected by requir- 
ing the highest momentum track in the event (excluding 
tracks from the -Btag reconstruction) to have a momen- 
tum of 1.7GeV/c < p* < 3.0GeV/c and to satisfy parti- 
cle identification (PID) criteria for either an electron or 
a muon. In events with a charged -Btag, the charge of the 
track is required to be opposite that of the -Btag, while 
for a neutral -Btag the high-p* lepton is permitted to have 
either positive or negative charge. 

Once the -Btag and the signal lepton candidate are iden- 
tified, -B+ I'^v events should ideally have no other par- 
ticles in the detector, while B^ — > t^T~ events should 
additionally contain only the t~ decay daughters. For 
the latter, the t~ rest frame is calculated from the ob- 
served signal lepton, assuming the nominal energy and 
momentum of the for a 2-body -B" decay. The six t 
decay modes considered are listed in Table HI The sec- 
ond highest momentum track in the event (again, exclud- 
ing -Btag reconstruction) is assumed to be a r daughter, 
and is required to have a charge opposite to the primary 
signal lepton. If this track satisfies electron or muon 
PID, the event is considered to be a leptonic r decay. 
Otherwise, the track is assumed to be a pion and the 
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quantity A£^^ is calculated for the hadronic decay modes 
listed in Table HI AEr — J2^TT-y> + Pu — rrir, where 
rrir = 1.777 GeV/c^, the sum is over the r daughter candi- 
dates, the momentum of the neutrino is pi, = | -Pjt- .-n-o I > 
and all quantities are measured in the t~ rest frame. We 
assign the decay mode for which \AEr \ is smallest, requir- 
ing additional conditions for the decay modes that pro- 
ceed through the intermediate resonances p~ —> tt^tt*^, 
ai~ —^ TT^TT^Tr" and ai~ ■k~'k~tt~^. We calculate 
the quantity cos6't— p — {2ErEp — — mp)/(2|p^||p^|), 
where {Et^Pt) and {Ep,pp) are the four-momenta in the 
^signal frame, and mr and mp are the masses of the 
T and p. For a correctly reconstructed p, this quan- 
tity peaks near unity. If the candidate does not satisfy 
cos 9r-p > 0.70 the mode with the next smallest \AEr \ (if 
one is present) is selected instead. Analogous quantities 
are calculated for the t~ — > tt^ tt*' tt^ Vt and r~ — > tt" 
7r~ TT"*" Vt modes, but with an instead of a . The 
requirements of cos 6*^-01 > 0.45 and cos^t- ai > 0.35 
are used for the two cases, respectively. There are no 
additional requirements on the p or ai. 



momenta. The signal events are selected by requiring 
APmiss to be less than 0.5GeV/c. 

For t^v modes we also consider the direc- 

tion of the missing momentum cos9p^_^:^^^ = Pzmiss/Pmiss, 
where the subscript z indicates the component of the 
momentum in the direction parallel to the beam pipe, 
as measured in the T{AS) CM frame. The requirement 
—0.76 < cosOp^.^^ < 0.92 is determined by the geometry 
of the detector; events where Pmiss points outside of the 
detector acceptance in the forward or backward direction 
are excluded. 

The quantity E'cxtra = J2 -E^track + J2 -E^clustor - Ei,+ - 

J2 -^e- ,7T- ,ttO describes the amount of energy recorded by 
the detector that is not accounted for by the high mo- 
mentum lepton and t~ daughters (in the case of B'^ — > 
T^). The clusters and tracks associated with the re- 
construction of Stag are excluded from the sums, and 
only clusters with energy more than 50MeV in the CM 
frame are considered. We require i?extra to be less than 
1.0 GeV in the CM frame. The signal and background 
distributions for -Ecxtia are shown in Fig. [TJ 



TABLE I: The r decays considered are listed with their 
branching fractions, in percent [l8| . 



r decay mode 


Branching Fraction 




17.84±0.05 




17.36±0.05 


7r~ fr 


10.90±0.07 


TV~ 7r" Ut 


25.50±0.10 




9.25±0.12 


TX~ 7r~ TT^ Ur 


9.33±0.08 



Additional background, for both i?"*" I'^v and B^ 
t^T^ decays, can arise from a variety of sources, includ- 
ing beam backgrounds, unassociated hadronic shower 
fragments, reconstruction artifacts, bremsstrahlung and 
photon conversions. We demand that events have no 
more than two extra charged tracks and six extra neu- 
tral clusters, allowing the presence of low energy particles 
not necessarily associated with the decay of the T{AS). 
Requirements on the missing momentum and extra en- 
ergy in the event are utilized to ensure that such particles 
are unimportant for the analysis. Since many of the fol- 
lowing requirements are optimized for each signal mode 
individually, we quote the approximate values. 

The extra momentum in the event is represented by 
APmiss = |pmiss + I]P«,7rL whcrc p£,7r are the momenta of 
the lepton or pion candidate(s) assumed to be recoiling 
against the neutrinos. The missing momentum is calcu- 
lated according to p,„iss = Pr(45) - - Pali, where 
Pall is the momentum of all tracks and clusters left after 
the i?tag reconstruction. APmiss is calculated in the rest 
frame of the parent of the neutrino(s), so that the miss- 
ing momentum balances the sum of other signal particles' 
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FIG. 1: Ecxtia distributions for the background simulation 
and data (left) and the signal (right) after most of the se- 
lection criteria have been applied. The upper plots are for 
l^u modes and the lower plots are for _B° — > 
modes. The background distributions show electron and 
muon modes together, as they are nearly identical. The back- 
ground is almost completely dominated by BB events. The 
signal modes are shown with a branching fraction of 10~^. 



The signal yields are extracted from unbinned maxi- 
mum likelihood fits to the signal lepton momentum distri- 
butions, as measured in the Psignai frame. The signal and 
background MC distributions are fitted by phenomeno- 
logical probability density functions (PDF). The si gnal 
distributions are modeled with Crystal Ball functions 21 1 
to account for the energy loss due to unreconstructed 
bremsstrahlung photons. The B'^ ^ £~^u background is 
modeled with an exponential decay and a Gaussian dis- 
tribution, while the B'^ £^ t" background is modeled 
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with a double Gaussian distribution. The PDF param- 
eters are determined from simulated events. The fit is 
performed using the following likelihood function: 



N 



nK/,(z) + nfc/fc(i)], (2) 



i=l 



where TV is the total number of events in the fit region, 
fs{i) and /&(«) are the PDFs for the signal and back- 
ground, and Ub and Hs are the number of background 
and signal events. All parameters of the signal and back- 
ground PDFs remain fixed, while Us and Ub are allowed 
to float. The fits are restricted to the ranges in p* shown 
in FigH 

The 90% confidence level (C.L.) upper limit on the 
branching fraction B is determined by solving for 6^0% 

in 0.90 = /g C{B)dB/ C{B)dB for events lying 
in the signal regions of 2.40GeV/c < p* < 2.75GeV/c 
for B+ £+iy and 2.20GeV/c < p* < 2.42GeV/c for 
— > £+T~ . B is related to the signal yield ti* through 
a substitution n* — etot x 2 x iVgg x B, where ctot is the 
total signal selection efficiency and -/Vgg is the number 
of B'^B" or B^B^ pairs in the data sample. The sig- 
nal selection efficiencies, expected number of background 
events and fit results are given in Table |TT1 The num- 
ber of signal events given by the fits is consistent with 
zero for all decay modes. The uncertainties in Table Hi] 
are statistical except for those shown for B which are the 
statistical and systematic uncertainties added in quadra- 
ture. 
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FIG. 2: The unbinned maximum likelihood fits on the lepton 
momentum. The dashed line, representing the signal PDF 
with an arbitrary scaling, indicates where the signal is ex- 
pected. 



The systematic uncertainties arising from the fitting 
procedure are studied by repeating the procedure on ad- 
ditional simulated samples, generated according to the 
PDFs, with varying number of signal events. System- 
atic effects are studied by repeating the procedure with 



TABLE II: Signal selection efficiency etot determined from 
MC, the fitted numbers of signal and background events in the 
signal regions and nl , and the branching fractions B. The 
uncertainties for B include statistical and systematic terms. 
The uncertainties for the other quantities are statistical only. 







Signal Mode 












etot X lO"" 


135 ±4 


120 ± 4 32 ± 2 


27 ±2 


nl MC 


2.66 ±0.13 


5.74 ±0.25 8.69 ±0.27 


12.14 ±0.45 


nl 


2.67 ±0.19 


5.67 ±0.34 9.35 ±0.35 


13.03 ±0.31 


n*s 


-0.07 ±0.03 


-0.11 ±0.05 0.02 ±0.01 


0.01 ±0.01 


B X 10-^ 


-o.it?:? 






g90% V.L. 


5.2 X lO"'' 


5.6 X 10"'* 2.8 X 10"^ 


2.2 X 10"^ 



PDF parameters varied by their uncertainties. For the 
case of zero signal events, we find negligible effects on the 
branching fraction values, and take the standard devia- 
tion of Us and Ub from their expected values in the fits 
as systematic uncertainties. We find the fits to be well 
behaved and having no significant sources of bias, intro- 
ducing no additional uncertainties. Total uncertainties 
associated with the fitting procedure are listed in Table 
III for each decay mode. 

The discrepancies between simulation and data are 
treated as detailed in the following paragraphs. The 
number of correctly reconstructed i?tag events in the toes 
signal region is compared between simulation and data. 
The Toes distributions for simulation and data are fit- 
ted with a combination of ARGUS [IJ] and Crystal Ball 
functions, allowing the number of TOes peaking events to 
be estimated by integrating the peaking component be- 
tween 5.270 GeV/c2 and 5.288 GeV/c^. We find the sim- 
ulation to underestimate the number of events with a 
good Stag and scale the signal selection efficiency by a 
factor of 1.11±0.06 (1.05±0.06) for events with a neutral 
(charged) Stag. 

In addition, the PID efficiencies in simulation are cor- 
rected for the 2-5% lower efficiencies in data. We assign 
associated uncertainties of about 2% for high momentum 
particles (signal lepton), and about 5% for tau daughters. 
The misidentification rate of leptons and pions is found 
to be negligible in the simulated samples, after all selec- 
tion criteria are applied. An uncertainty in the track- 
ing algorithm introduces an additional 0.8% systematic 
uncertainty for each charged track present in any given 
signal mode (e.g. 1.6% for B^ £^ , tt" Vr). 

The uncertainties for B'^ modes are calculated 

as weighted averages of all decay modes. 

Table IIIII lists the sources and the magnitudes of the 
uncertainties with their effect on B. The uncertainties are 
incorporated into the final results by varying the branch- 
ing fraction assumption by its uncertainty when integrat- 
ing C for the 90% C.L. upper limit. 
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TABLE III: The sources and magnitudes of systematic uncer- 
tainties, in percent. 
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Signal Mode 




Uncertainty source 


e+ r 


fX^ T 


e^v 




Signal Fit 


5.6 


10.6 


4.3 


8.2 


Background Fit 


3.9 


3.1 


5.1 


7.8 


Stag efficiency 


6.4 


6.4 


5.8 


5.8 


FID efficiency 


5.3 


5.8 


1.0 


2.0 


MC Statistics 


8.6 


7.4 


3.0 


2.8 


Tracking efficiency 


1.7 


1.7 


0.8 


0.8 




l.f 


1.1 


1.1 


1.1 



We have presented searches for the rare leptonic decays 
5+ C'^v and B'^ £^t^, where £ = e,fi, using a 
novel hadronic tag reconstruction technique. We find no 
evidence of signal in any of the decay modes in a data 
sample of approximately 378 million BB pairs (342 fb^^), 
and set the branching fraction upper hmits at B{B+ 
e+iy) < 5.2 x IQ-^, B{B+ fi+i/) < 5.6 x 10-^ B{B° 
e+T-) < 2.8 X 10-^ and B{B" M+t") < 2.2 x 10~^ 
at 90% confidence level. While these upper limits on 
B{B'^ — > e^iy') and B{B^ M^^) complement the more 
stringent limits available from inclusive studies 13, 1^ . 
the B^ e~^T~ and /i+r^ results are the most 

stringent published upper limits available. 
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